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Abstract 

Light experiences a non-uniformly moving medium as an effec- 
tive gravitational field, endowed with an effective metric tensor g'^'^ = 
^t^u _|_ ^j^2 _ i^^M^*^^ ji being the refractive index and the four- 
velocity of the medium. Leonhardt and Piwnicki [Phys. Rev. A 60, 
4301 (1999)] argued that a flowing dielectric fluid of this kind can be 
used to generate an 'optical black hole'. In the Leonhardt-Piwnicki 
model, only a vortex flow was considered. It was later pointed out 
by Visser [Phys. Rev. Lett. 85, 5252 (2000)] that in order to form 
a proper optical black hole containing an event horizon, it becomes 
necessary to add an inward radial velocity component to the vortex 
flow. In the present paper we undertake this task: we consider a full 
spiral flow, consisting of a vortex component plus a radially infalling 
component. Light propagates in such a dielectric medium in a way 
similar to that occurring around a rotating black hole. We calculate, 
and show graphically, the effective potential versus the radial distance 
from the vortex singularity, and show that the spiral flow can always 
capture light in both a positive, and a negative, inverse impact pa- 
rameter interval. The existence of a genuine event horizon is found 
to depend on the strength of the radial flow, relative to the strength 
of the azimuthal flow. A limitation of our fluid model is that it is 
nondispersive. 
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1 Introduction 



The classical topic of light propagation in moving dielectric media has re- 
ceived considerable interest in recent years - for a review of recent develop- 
ments see, for instance, Ref. Much of the reason for the revitalization 
of this field of research is evidently the advent of the electromagnetically 
induced transparency (EIT) technique 0, whereby the light velocity in mat- 
ter can be slowed down dramatically, to orders of 1-10 m/s |@. Under such 
circumstances it becomes possible to generate a fluid flow whose velocity is 
comparable to that of light in the fluid. The flow will then have an important 
influence upon the propagation of light. 

Light waves experience a moving dielectric medium like an effective grav- 
itational field. The effective metric, g'^" = rj^^ + (n^ — l)u^u^ , was used 
long time ago to develop a relativistic model of a non-dispersive dielectric 
medium, and will be applied also in this paper, n is the refractive index, and 

is the four-velocity of the medium. The use of the effective metric allows 
us to apply the general-relativistic formalism to systems that normally are 
not identified with general relativity. Some earlier papers along these lines 
are Refs.[4-7]. Recent discussions on the general-relativistic analog models 
can be found at the web-address of the Rio de Janeiro Workshop in October 
2000 p. Since the velocity of light can be slowed down so dramatically, 
the flowing medium gives rise to the effect that light rays are moving along 
curved trajectories. 

Central references in our present context are the papers of Leonhardt 



and Piwnicki [0, ITO[. These authors considered a dielectric fluid swirling 



around a straight linear vortex, arguing that such a system can be used to 
generate an 'optical black hole'. Subsequent remarks of Visser made it 
clear, however, that the interpretation of the Leonhardt-Piwnicki model in 
its original form as an optical black hole is not feasible: in order to obtain a 
black hole endowed with an event horizon one has to add an inward radial 
component to the vortex flow (cf. also the rebuttal |]T2[|). This brings us to the 
subject of the present paper: We will develop the optical black hole theory 
based upon a general spiral flow model, consisting of a vortex component 
plus a radially infalling component. We introduce an "effective potential", 
constructed as the appropriate generalization of the conventional effective 



potential outside a Schwarzschild black hole |T^. Plotting the variation of 
the effective potential versus radial distance r from the vortex singularity, we 
get a powerful tool for predicting particle and light trajectories. Our explicit 
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calculations prove that Visser's remarks |]TT[] (which are based upon general 
considerations) are right. We find that the spiral fiow can always capture light 
within a positive, and a negative, interval for the inverse impact parameter. 
If the azimuthal velocity component is increased, the parameter interval for 
light being able to escape from the system is also increased. This reminds 
us of gravitational black holes in that way that it is known that a rotating 
black hole captures light less efficiently than a static one of same mass. 

One important limitation of our theory should be emphasized: we con- 
sider non- dispersive media only. That means, we are unable to cope with the 
extremely dispersion-sensitive EIT effect [^. Some more comments should 
be given, concerning the physical background for this statement. Although 
the EIT technique is extremely effective - the quantum interference allows 
pulse transmission through an atomic cloud which would otherwise have a 
transmission coefficient only of order e~^^° - it is notable that the effect is 
also extremely dispersive. The steep dispersive profile is shown by Fig. 2 
in the article of Hau et al. According to the figure the high refractive 
index prevails only over a frequency interval of about Az/ = 2 MHz around 
the central frequency v = 0.5 x 10^^ Hz (corresponding to a wavelength of 
A = 589 nm). When the light propagates in the direction of motion of the 
fiuid, it follows from a consideration of the four-momentum of the pho- 
ton that the fiuid velocity v gives rise to a Doppler shift of Au = unv/c, 
OT V = (c/n)(Az//z/) = 120/72 cm/s. Now, in the atomic vapour n is prac- 
tically unity, so we see that already a moderate fiuid velocity of about 100 
cm/s is sufficient for the Doppler effect to bring the frequency in the local 
rest inertial frame of the fluid outside the EIT-generated high-?7, band. The 
Doppler argument assumes for simplicity that n is a constant, but this is not 
disturbing for our estimate: both is n in the undisturbed vapour close to 
unity, as mentioned, and moreover the refractive index-disturbance caused 
by the laser beams is very small. An = 1.8 x 10~^. In practice, dispersive 
effects thus play an important role; the Doppler detuning will exceed the EIT 
window, and the present theory is not directly applicable to slow light. We 
think, however, that the nondispersive theory is of interest in its own right. 

In the next section we summarize the basic properties of the phenomeno- 
logical electromagnetic theory in a dielectric medium, emphasizing the formal 
analogy to general relativity, and also the connection between this kind of 
theory and the transformation procedure worked out some years ago 0, 
from which one can map the electromagnetic fleld in a vacuum to that in 
a medium. Then, in Sec. 3 we consider the effective metric for the two 
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components of the spiral flow. In Sec. 4, before embarking on the rather 
complex theory of the general spiral case, we define the "effective potential", 
which is a new concept in this context, and apply it to the two basic flow 
components separately (the pure sink, and the pure vortex). Figures 2 and 
3 show the "effective potentials" for these two cases. The figures are useful, 
in that they show which regions are inaccessible for photon trajectories. In 
Sec. 5, which is the main section, we work out the "effective potential" for 
the full spiral flow, and introduce various concepts: the event horizon, the 
apparent horizon, the ergo region, and the static limit. Various kinds of pho- 
ton trajectories are discussed in relation to Fig. 4. As far as we know, this 
kind of analysis has not been given before, for the general composite flow. A 
summary is given in Sec. 6, along with a closer discussion on the relationship 
to a Schwarzschild black hole. 

We put c = 1, adopt the Minkowski metric convention rj^^^ = (1,-1,— 1,-1), 
and follow the basic formulation of Leonhardt and Piwnicki |^. We men- 
tion in passing that some other related recent treatises can be found in Refs. 
[14-19]. 

2 Basic formalism 

We begin by giving a brief account of the basic formalism, letting in this 
section the four-velocity u'^ of the fluid be unspecified (thus not necessarily 
being the sum of vertex flow and radial flow). With the above metric, u^u^ = 
1. If W are the components of the three-dimensional velocity vector U, i.e., 
U^Ui = U^, we have thus 

«^ = 7(l,t/^ (1) 

where 7 is the Lorentz factor, 7 = (1 — U^)"^''^. 

It ought to be emphasized that we are applying a kind of hybrid interpre- 
tation to the system. First, as far as the fluid flow velocity m'^ and the wave 
four-vector are concerned, we consider the system as an ordinary special- 
relativistic system, in which covariant and contravariant indices are moved 
by means of the Minkowski metric t]^^. Second, when considering the wave 
propagation through the moving flow, we make use of the effective metric, 
which will be designated by ^^j,. The background for introducing this kind of 
metric can be seen as follows. Starting from Maxwell's equations for a pure 
radiation field in a nondispersive medium moving with special-relativistic ve- 
locity we obtain, by defining the four-potential from the basic field 
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tensor F^^, via the relationship F^^, = d^A^ — d^A^, the field equation |^ 

[□ + K{d.uf]A^ = (2) 

(□ = dfj,d'^, d.u = d^u^), where we have defined k as 

K = n^ -I. (3) 

In momentum space, Eq. (2) implies 

+ K{k ■ m)2 = (4) 

(A;^ = k^k^). Since the indices of the photon four-momentum fc^, as men- 
tioned, can be lowered and raised by the Minkowski metric t]^^^, we can write 
Eq. (4) as 

g'-'k^K = 0, (5) 
with gf^" being the effective metric 

gf^-^ = rit^- + (6) 

It is the appearance of Eq. (5) that makes it possible to borrow terminology 
from general relativity and consider the propagation of light waves as taking 
place in a fictitious or 'effective' curved space whose metric is given by Eq. (6). 

In the limiting case of Hamilton's geometrical optics, which is adopted 
here, the electromagnetic field tensor is modelled by the eikonal ansatz ||20l . 



It implies that the photon four-momentum is essentially the four- divergence 
of the eikonal S, i.e., /c^ = — 9^5*. Accordingly, it becomes natural to write 
the wave equation in such a way that it involves the covariant components 
of the wave vector, as in Eq. (5). 

We mention in passing the following point. In Ref. we developed 
a transformation procedure which maps the electrodynamic theory for a 
vacuum field directly into the phenomenological nondispersive theory for 
a medium field. This mapping procedure is actually quite powerful, espe- 
cially in connection with quantization of the medium field. Consistent with 
the use of the canonical (space-like) four-momentum fc^ in the medium, the 
quantization is done in accordance with the nonsymmetric, divergence-free 
Minkowski energy-momentum tensor. The central building block in this kind 
of theory is the matrix b'^'^ , defined by 
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satisfying the relationship (p an integer): 

= n^''' + (nP - (8) 

The effective metric in Eq. (6) is thus nothing else than the second power 
of the basic matrix: g^" = {b'^Y^ . 



3 The two flow components 

We now turn to the situation in which the flow, as mentioned above, is com- 
posed of a vortex component and a radial component. The three-dimensional 
fluid velocity U can thus be written as 

U = U'e, + f/^r)e^, (9) 

in an orthogonal basis (designated by hats) {e^} = {ef,e^,ez}. We write 

the azimuthal component in the form f/*^ = Vt/r with Vt a constant, r being 
the radial distance from the singular vortex (situated along the z axis). In a 
coordinate basis, {ej} = {d/dr, d/dcj), d/dz}, this corresponds to = fl/r"^. 
The relativistic correction is introduced via the Lorentz factor 7, so that 
If^ — > f2/(7r^). As for the radial flow, we start from the nonrelativistic 
standard form for a singular line source, = T/r, T being a constant that 
is negative when, as assumed here, the source is actually a sink. (There is 
for the radial component no difference between the orthogonal basis and the 
coordinate basis, i.e., = U"^.) Again including the relativistic correction, 
we have —>■ r/(7r). Altogether, the three-dimensional fluid velocity U in 
coordinate basis becomes decomposed as 

u = ([/^f/^f/^) = (-,-^,o), (10) 

with Lorentz factor 



i = \li + ^^. (11) 



From Eq. (1) we obtain the fluid's four-velocity: 

«^ = («^«^«^«^) = (^7,^,^,0). (12) 
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There are two nontrivial Killing vectors in the problem, viz. C,{t) = d/dt and 
^(0) = d/d(j). The four-vector products k ■ ^(t) and k ■ ^(0) are constant along 
the photon trajectory, meaning that the components ko and fc^ are constants. 
We will designate them respectively by E and —L: 

ko = E, k^ = -L. (13) 

The physical meaning of these quantities can be seen as follows. Consider 
the line element for light at spatial infinity: 

ds"^ = g^^dxf'dx'' ^dt^ - dr^ - r'^d(j? - dz^ . (14) 

In an orthonormal basis {uJ^} at r — oo we thus have that the zeroth 
component becomes uj^ = dt/n. The local photon energy Eio^ai at infinity is 
the contraction of the four-momentum k^ with a;°: 

- ~ 

Eiocai = k) = — = nko, (15) 
n 

since g^^ = rt^. Thus 

E=^l^. (16) 

n 

Since E and Eiocai differ only through a constant factor n, we can use these 
two concepts interchangeably. 

As for the azimuthal part, the interpretation of L can be taken over di- 



rectly from the theory of an ordinary black hole [|I^ . The spatial components 
of the metric at infinity are identical to those of the Minkowski metric, for 
the Schwarzschild metric as well as for the present case with effective met- 
ric g^'^ . The quantity L is the angular momentum (for an ordinary black 
hole referring to equatorial orbits). The minus sign in front of the second of 
Eqs. (13) is due to our metric convention. 
The quantity 6, defined by 

will henceforth be referred to as the impact parameter (at infinity). (This is 
not in perfect agreement with the usual definition of the impact parameter.) 
As discussed by Leonhardt and Piwnicki |^ the canonical angular momentum 
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L differs from the kinetic angular momentum by a constant, due to the optical 
Aharonov-Bohm effect. 

The light paths are defined by the wave four-vector tracing out zero 



geodesies of g'^'^. Introducing, as usual |13], an affine parameter A, we can 



parametrize the light path through the relation 



Of main interest for us is the r— component: 

It is to be noted that when dr/d\ is zero the photon, with its radial mo- 
mentum equal to k"^ , has a turning point. 

For reference purposes we write down the components of the effective 
metric for the full spiral flow. From Eqs. (6), (10), and (11) we have 

00 - ' - " -t- i 



= 1 + K 1 



-Or- = ^ /l + + 



r = -1 + — , (20) 



y ^3 ) 



1 kVI 
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r = -1- 

We will need also the covariant components: 

fl-oo = 1 ^ 1 + 



^2 I J.2 
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90r 



9o^ 



Tn?r 



1 + : 



nil I n^ + T^ 

1 ~\~ I^^ • 



(21) 
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_^2_ 



/tfi2 



-1. 



Op. 



These expressions satisfy the contraction conditions g^'^g^j.p 

It ought to be pointed out here that we have constructed the velocity of 
the composite flow, Eqs. (10)-(12), as an ad hoc generalization of the flow 
adopted by Leonhardt and Piwnicki 0, |1^. Since the four- velocity has 



to satisfy the relativistic Euler equation, the value of the pressure p has to 
adjust itself accordingly. We do not consider the pressure here. It should also 
be mentioned that the values of the metric components given by Eqs. (20) 
and (21) are different from those used by Linet |TB|. 



4 "Effective potential" for the spiral flow. Ap- 
plication to the special cases 

The Hamilton- Jacobi equation, Eq. (5), will now be solved for the general 
planar centrally symmetric flow. This spiral flow model of course gives asso- 
ciations to a Schwarzschild black hole. We will show that the system indeed 
gives rise to effects similar to those occurring near rotating black holes. Thus 
the concepts of an event horizon, and an ergosphere, can be introduced in a 
natural way, and we will also show how qualitative features of light trajec- 
tories in the flow are obtained by plotting the "effective potential" for light. 
First, the expression for the "effective potential" has to be derived. 
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4.1 The "effective potential" for the spiral flow 



Employing the effective metric picture, we will use g^" to raise indices of the 
covariant wave ve 
will be defined as 



covariant wave vector A;^. Thus the contravariant effective wave vector 



= gt'-^k^^ (22) 

the inverse equation being /c^ = g^vk'^ . (Note that the covariant vector A;^ 
retains its usual special-rclativistic meaning.) In order to solve Eq. (5), one 
needs to find all the components of expressed in terms of known quantities. 
Since three components of k^ are constants, /cq — k^ — —L, k^ — 0, it 
remains to find the component kr- To this end we first express kr as 

kr = (jrok^ + Qrrk'' + gr<i>k'^ ■, (23) 

and thereafter insert k^ = g^'^k^, k'^ = g'^'^k^ to convert Eq. (23) to an 
expression in which all components except from kr are known. Taking into 
account the constants above, we can then write kr as 

■A ^ B ^ Orr dr , V 

where 

B ^ ~gor~g'^ + ~g^r~9'^t (25) 



C=l~~gorr-~9^r~g'^''- 

At the turning point of a trajectory dr/dX = 0, so that the last term in 
Eq. (24) vanishes. Turning points will be a matter of main interest below. It 
is convenient, therefore, to introduce a new symbol kr denoting the value of 
kr at just the turning point. It is given by 

K^^E- |l. (26) 

All components of k^ in the Hamilton-Jacobi equation are thereby known. If 
we for simplicity consider the solution of this equation accepting the condition 
(26), the Hamilton-Jacobi equation determines which turning-point value of r 
is compatible with given input values of E and L. To emphasize notationally 
this particular case, we give an extra 'hat' to the symbols, so that E 
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E, L ^ L. Whereas for given values of E and L the equation determines 
for which r the turning of the photon trajectory takes place, it should be 
noted that at the same r other photons, corresponding to different values of 
E and L are also allowed to exist, but they do not have turning points at 
this position. 

As will be discussed in more detail below, the turning point solutions 
give the effective potential for the system. The Hamilton- Jacobi equation is 
a quadratic equation in E (or alternatively in L); it is convenient to work 
instead in terms of the ratio between these quantities, i.e., = E/L, which 
is the inverse impact parameter. We write the Hamilton- Jacobi equation in 
the form 

p + Q + 7^ = o, (27) 

with known coefficients V, Q, TZ. This equation has two sets of solutions, 
61,^^ (r) and bZ^{r), which in the following will be called the "effective poten- 
tials" . 

Some remarks are called for here, since we are using this concept in a 
somewhat more restricted sense than what is usually the case. Recall the 
conventional equation determining the azimuthal angle (p: 

{L/r'^)dr 



v/2m[Eo-K//(r)]' 

m and Eq being constants. In this case the r-dependent function Ve//(r) is 
naturally called an effective potential. For example, trajectories around a 
Schwarzschild black holes can be expressed in this form. The physical mean- 
ing of our "effective potential" is simply that it corresponds to a curve from 
which the turning points for photons can be read off. In other words, whereas 



in Ref. |T^, for instance, the effective potential is part of the equation of 
motion, regardless of whether the trajectory has a turning point or not, in 
our model the "effective potential" is found by solving the equation of mo- 
tion precisely at a turning point. Our "potential" moreover does not have a 
nonrelativistic analogy. 

The expressions for V, Q, TZ are at first sight complicated, but after a 
lengthy algebra making use of the metric components (20), (21) as well as 
Eqs. (25), (26) we arrive at the following relatively simple expressions: 



V 
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2-^^V^ + ^> (28) 

The solutions of Eq. (27) now become, when we make use of the expressions 
(28), 



4.2 Inaccessible regions 

In order to justify our interpretation of the expression (29) as an "effective 
potential", we shall now show how it acts an an impenetrable barrier. Let 
us assume temporarily that dr/d\ is nonzero. Then kr, instead of kr given 
by Eq. (26), has to be inserted in Eq. (5). The Hamilton- Jacobi equation 
accordingly becomes 



r%?(^V, ■ (30) 



C d\ 



with the same factors (25) and (28) as before. Here has taken the place 
of h^^. since turning points solutions arc now not searched for. 

Consider the left hand side of Eq. (30). From Eq. (28) it follows that if 

> «;F^ (31) 

then P > 0. If the left hand side is plotted versus the curve gets a positive 
curvature, and crosses the abscissa axis at the turning point solutions, 6^^, 
as indicated in Fig. 1. 

Consider next the right hand side of Eq. (30). The first term turns out 
to be zero, as it must be, since the right hand side can not depend on the 
sign of dr/dX while the left does not. Thus, Eq. (30) becomes 

!,-2 , /-I L-l , i:. ~rr{9rrY f dr 



vh-' + Qh-' + n = -r'-^[^] . (32) 
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From Eq. (20) it follows that g"^"^ is negative in the outer region defined by 
Eq. (31), so that the right hand side of Eq. (32) is positive in this region. The 
condition that both sides of Eq. (32) must be positive, makes us thus conclude 
that the region between bZ^ and b^^ in Fig. 1 is inaccessible. The solution of 
the Hamilton-Jacobi equation (32) accordingly contains a potential barrier. 



Left 
side 








Figure 1: Qualitative sketch of the left hand side of Eq.(32) versus b ^, when 
> kF^ 

Let us now analyse the inner region, 

< kT^. (33) 

From Eq. (29) it follows that there are no real solutions for turning points 
in this region; the left hand side of Eq. (30) has no zeros. The inequality 
(33) implies that V is negative, so that the curve in Fig. 1 is to be replaced 
by a curve (not shown here) which is concave downwards and which has no 
intersection point with the abscissa axis. Since g^^ is now positive, the left 
hand side of Eq. (30) is negative for any value of the impact parameter. Any 
photon, if it comes that far, is allowed to exist for any value of r in the 
inner region. There is no potential barrier, so that a photon approaching the 
singularity will ultimately get absorbed there. 

4.3 Application to the perfect sink flow 

Before embarking on an analysis of the full spiral flow, it is worthwhile to 
consider the two special cases flrst. 
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The simplest special case is that of a perfect sink, corresponding to f2 = 
0, r < 0. According to Eq. (29) the effective potential is now 



(34) 



The "effective potential" is shown qualitatively in Fig. 2. The symmetry 
of the curve about the abscissa axis reflects the physical property that the 
turning point of a photon with angular momentum L must occur at the same 
value of r as the turning point of a photon with angular momentum —L. 




Figure 2: Perfect sink: effective potential, = E/L, versus r. Broken lines, 
i?2, Rs and i?4 refer to the three possible photon-trajectories. 

The present "effective potential" is essentially the same as the effective 
potential for light in Schwarzschild geometry [|r^ . At the position determined 

by 

= kF^ (35) 

the effective potential becomes zero. This point is interpreted as the event 
horizon of the system. At smaller r there is not, for any impact parameter, 
a turning point. 
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The fluid's three-dimensional velocity in coordinate basis is, according to 
Eq. (10), 

U = (^^, 0, 0^ , with 7 = ^1 + ^^. (36) 

Now, the event horizon would from a physical viewpoint be naturally deflned 
as the position where the velocity of light relative to the medium, 1/n, is 
equal to the inward radial fluid velocity, T/'jr. This leads to the condition 
l/n"^ = r^/7^r^ which, upon insertion of the Lorentz factor from Eq. (36), is 
just the same as Eq. (35). This shows the consistency of identifying the zero 
point of the effective potential with the event horizon . 

The effective potential curve for a perfect sink is essentially the same as 
for light outside a Schwarzschild black hole. Whereas a conventional black 
hole can be visualized as a sink of space-time, the present 'optical black hole' 
can analogously be visualized as a sink of the moving medium, which at the 
same time serves as the artiflcial 'space-time'. 

The three horizontal broken lines marked by R2-,Rz-,Ri in Fig. 2 refer 
to the fact that there are three different types of motion for photons in the 
radial flow. It is here flrst of all appropriate to recall that in Schwarzschild 
geometry there are actually four types of motion possible when the orbiting 
particle has nonzero mass The flrst type, which we would naturally 

denote by refer to the 'pit' in the effective potential that permits stable 
orbits to exist around the black hole. Now, there are no stable orbits for 
light, so that Region Ri does not exist, neither in Schwarzschild geometry 
nor in 'perfect sink' geometry. 

The other three regions exist for light, in Schwarzscild geometry as well 
as in the present case: 

Light can approach the centre of the optical black hole with an inverse 
impact parameter that lies above the "effective potential". If h^^ at some 
radius is equal to the "effective potential" 6^^, this will be the turning point 
for the trajectory. After turning, the light goes outwards and escapes to 
inflnity. This is illustrated as Region 2. 

If the photon has an inverse impact parameter that lies above the max- 
imum of the "effective potential", it will be captured by the optical black 
hole. This is shown as Region 3. The region also covers the case where a 
photon, emitted somewhere outside the event horizon, escapes to inflnity. 

Region 4 is the region for a photon which is emitted somewhere near 
the event horizon and eventually is lost to the inside region. It may at flrst 
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propagate in the outward direction, but 'hits' the potential curve to the left 
of the maximum point, and turns back into the centre of the optical black 
hole. 

Finally, we comment briefly on the unstable photon orbits that may exist 
temporarily in the pure sink flow. Differentiating Eq. (34) we see that the 
extrema of the "effective potential" are lying at 

= 2kT^. (37) 

The two extrema are lying at the same position r (which is physically obvious 
since they describe photons travelling in opposite directions past the opti- 
cal black hole), and they correspond to unstable orbits. The corresponding 
extrema follow from Eq. (34) as 

= ± (38) 

Evidently such orbits will never be permanent; a photon having impact 
parameter very close to the critical one given by Eq. (38) will spiral several 
times around the centre before it either escapes to infinity or else becomes 
lost inside the event horizon. 



4.4 Application to the perfect vortex flow 

This is the case studied by Leonhardt and Piwnicki 1^, [l^. It corresponds 
to r = 0. We will assume a positive vortex, so that f2 > 0. From Eq. (29) 
we get the "effective potential": 

^ t,n^l + n^/r^±nr ^^^^ 

A qualitative plot of the effective potential is shown in Fig. 3. The rotation- 
induced asymmetry is obvious. The upper curve corresponds to the positive 
sign in Eq. (39). When r approaches zero, both solutions diverge. As in the 
general case, the region between the two solutions is inaccessible. 

When VL is positive, the motion is in the positive direction. Positive 
rotation of the fluid corresponds to positive angular momentum, so that 
photons with L > move along the medium whereas photons with L < are 
fighting against the flow. This is the physical reason why photon trajectories 
corresponding to L and —L are non-symmetric in Fig. 3. 
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Figure 3: Perfect vortex: effective potential, b^^ = E/L, versus r. Broken 
lines, R2, R3 and R4 refer to the three possible photon-trajectories. 

For later use, let us introduce the ergo region concept. It will be showed 
below that the boundary of this region is at the point where the negative 
solution bZ for the "effective potential" in Fig. 1 coincides with the origin 
of the b~^ axis: 

bz\r) = 0. 

From Eq. (39) this means 
As from Eq. (10) we have, for vortex flow, 

u= (0,^,0^ 

with 

(corresponding to the orthogonal azimuthal component = Q/{'jr)), we 
see that the same equation (41) can alternatively be found by solving the 
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(40) 
(41) 

(42) 

(43) 



relationship fi^/(7^r^) = Xjr?. This relationship means that the negative 
solution for the effective potential crosses the r axis in Fig. 3 where the 
(orthogonal) flow velocity f2/(7r) just counterbalances the velocity of light 
1/n relative to the medium. Inside the radius given by Eq. (41) the flow 
velocity exceeds the local velocity of light, and the region is per definition an 
ergo-region. The physical consequence of this is that light in this region is 
forced to participate in the positive 0-rotation around the centre, no matter 
which way the light tries to escape. 

The negative solution of Eq. (34) has a bottom point. This point corre- 
sponds to unstable circular orbits 0. One might imagine that light, having 
a value of just corresponding to a turning point at the bottom point, 
does not go in orbit at all, but stands still. If this scenario were true, light 
would go directly against the flow at a radial value where the flow velocity 
equals the velocity of light in the medium. 

This idea of still-standing light is however not correct. As was showed 
above, the point where the flow velocity equals the velocity of light is where 
the negative solution of the "effective potential" crosses the r axis. The 
bottom point of this potential is at a larger r— value, and the conclusion 
that we have to make, is that the bottom point represents an unsteady orbit 
instead. 

As pointed out by Visser [^, an azimuthal vortex cannot be taken to 
represent a genuine black hole. This can easily be seen from the "effective 
potential" plot of a photon possessing positive angular momentum. For, 
if the photon has a positive b~^, its accessible region is of type R2 (cf. also 
Fig. 2), meaning that the photon ultimately gets reflected out to infinity. The 
outward turning point can be found arbitrarily close to the vortex singularity. 
The system has no analogy to a Schwarzschild black hole in the sense of a 
point of no return existing for all photons. 

This does not mean that light cannot be captured by the vortex. As 
seen from Fig. 3 there is a region beneath the lower curve corresponding to 
type R3 (cf. again Fig. 2). In this photon coming in from infinity is 

captured by the vortex, and will never escape. This region exists only for 
light of negative angular momenta (and therefore of negative b~^). 

Photons having negative angular momenta, thus fighting against the 
stream, have no outer turning points inside the radius determined by the 
bottom point oi bZ^. The photons find the bottom point in some sense to be 
similar to an event horizon. To be a genuine horizon, however, the r— surface 
must have the property that no photon emitted on the inside has any possi- 
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bility to escape, even if emitted only slightly inside and in a radial direction. 
The bottom-point radius (or the radius in Eq. (41)) does not satisfy this 
condition, the reason being that the flow has no radial component. Photons 
can be captured by the vortex, since the flow at an early stage manipulates 
the trajectories and generates an inward directed spiraling motion for the 
photons. There is however nothing in the flow itself which drags photons 
directly towards the centre. If a photon in the inside region is emitted in 
a straight outward direction, it will continue to increase its r— value until it 
escapes to infinity. 

It is instructive in this context to recall the effective potential for Schwarzschild 
geometry: Photons coming in from infinity and pass the potential's top point, 
has no outer turning point and will be captured. This does not mean that the 
radius of the top point corresponds to the event horizon; the Schwarzschild 
radius lies at a lower r— value. 

Photons with negative can escape, if the absolute value of is less 
than the absolute value of the effective potential at the bottom point. Both 
types of regions, R2 and R^, exist for negative values of b~^, as shown in 
Fig. 3. 



5 Behaviour of light in the general flow 

Combining the two elementary flows considered above, we arrive at a con- 
siderably more complex physical picture. Figure 4 shows the effective poten- 
tial for the flow when the azimuthal and radial components are given equal 
weight: = T^. It is striking how the previous symmetric perfect-sink di- 
agram from Fig. 2 now becomes skew because of the addition of the vortex 
flow component. 



5.1 Event horizon 

The event horizon must first be properly defined. In general relativity it is 
defined as the boundary of the region from which null geodesies can never 
escape. If the inward-pointing radial component of the fluid velocity exceeds 
the velocity of light relative to the medium, then, no matter what direction 
the light ray tries to escape, it will be swept inward by the sink and trapped. 
This defines the event horizon as the boundary of superluminal radial flow. 
(An equivalent definition was made by Visser ||2^ for the event horizon in 
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Event Apparent ' static Limit 

Horizon Horizon (Outer limit for ergosphere) 



Figure 4: Effective potential, b'g^ = E/L versus r for light in a medium with 
both radial and azimuthal flow. Broken lines, R2, R3 and R4 refer to the 
three possible photon-trajectories. 



acoustics.) 

The mathematical definition of the event horizon is thus 



(cf. Eq. (10)). Taking into account Eq. (11) we get 

r = r event = VkT^ - (45) 

This equation shows that if the strength of the sink for a given vortex is so 
weak that 

|r| < (46) 

then Eq. (45) cannot be satisfied. It is noteworthy that the existence of 
an event horizon thus relies upon the relative strength of the sink/ vortex 
components, not upon their individual absolute magnitude. 

The event horizon is indicated by a vertical line in Fig. 4. A photon 
existing to the left of this line can never appear on the right hand side. 
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5.2 Apparent horizon, ergo region, and static limit 

As already mentioned, a pure sink acts as a static optical black hole and 
does not distinguish between light having positive or negative angular mo- 
menta. Assume now that the sink strength F is kept constant while the 
vortex strength Q gradually increases. The vortex nature of the flow pro- 
duces turning points 'higher up' for light having positive L while it 'takes 
away' turning points having negative L. Light propagating along with the 
azimuthal flow is helped to escape; light fighting against the flow is more 
easily captured. 

When Q is further increased the vortex nature gradually takes over, and 
a plot of the effective potential will be more and more like that of a pure 
vortex. There remains, however, the following significant difference. When 
there is a radial flow component, no matter how small in comparison with 
the azimuthal flow component, the surface where 

= (47) 

will have a special meaning. Inside this surface, which does not depend on 
Q at all, there are no turning points for incoming photons. This surface has 
much in common with an event horizon (it becomes just the same as the 
event horizon in the case of a perfect sink; cf. Eq. (35)), but is in general 
different from it. We will refer to the surface determined by the condition 
(47) as the apparent horizon. This definition of the apparent horizon stems 
from its visual apparance on the effective potential plot. It must not be 
confused with the definition made by Visser in of the apparent horizon 
for acoustic waves. 

The physical meaning of this concept is the following. Assume, as men- 
tioned, that r is constant while fl gradually increases. Plotting the "effective 
potential" versus r in the various cases one will see that the azimuthal flow 
influences the sink flow in much the same way as the rotation of an ordinary 
black hole influences the black hole's ability to capture light. It is known in 
general relativity that a rotating black hole captures light less efficiently than 



a static hole of the same mass |2^. The same kind of effect is found here. 
If r is kept fixed while fl is increased, the band of light that will escape 
increases. This band is the vertical distance between the top- and bottom- 
point of the effective potential. 

The vortex increases the band, but there is a limit to it. If ^ oo, the 
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last term in Eq. (29) becomes negligible and we get 



bg\fl ^oo)^P (48) 

the positive and negative solutions coincide. Here the meaning of the ap- 
parent horizon comes in: from Eq. (29) we see that the last term becomes 
imaginary outside r = rapp, where 



r, 



app 



v^|r| (49) 



is the radius of the apparent horizon. The "effective potential" is thus cut off 
at the apparent horizon. From Eqs. (48) and (49) we see that the maximum 
value of for a photon that can escape from the flow, no matter how large 
Q, is chosen, can be defined by 

i^app I r I 

With infinite f2, the lower curve for the "effective potential" does not cross 
the abscissa axis, so the maximum interval is < 6^^ < l/r^pp. 

The apparent horizon is marked by a vertical line in Fig. 4. Geometrically, 
the apparent horizon is at the r- value where the upper curve for the " effective 
potential" meets the lower, and the curves 'cut each other off'. The term 
'apparent' is used, because from a look at the potential plots, this radius 
appears to be the horizon. The apparent horizon is the point of no return 
for light coming in from the right on the figure. Inside this horizon there are 
no turning points, and light entering the region from the outside will never 
escape. Recall that the apparent horizon docs not depend on $7, in contrast 
to the event horizon, Eq. (45), which depends on both O and P. The event 
horizon always lies inside the apparent horizon. 

The difference between the apparent horizon and the event horizon is 
experienced by light that is emitted in the inside region. If this region extends 
beyond the apparent horizon but not beyond the event horizon, it is possible 
for the ray to escape to infinity. By contrast, in ordinary Schwarzschild 
geometry, and for the perfect sink, the light pulse must be emitted outside 
the apparent horizon if such an escape is to happen, since in these cases the 
apparent horizon coincides with the event horizon. 

The next concept that it becomes natural to introduce for this special 
kind of flow is the ergo region. In general relativity, a rotating black hole is 
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known to have an ergo region that reaches outside the event horizon. The 
part of the ergo region that hes outside the horizon is called the ergosphere. 
A characteristic property of the ergosphere is that qqq becomes negative, so 
that no state of rest can be defined. Correspondingly, the time Killing vector 
^(f) becomes spacelike, leading in turn to peculiar effects like the Penrose 



process (cf., for instance, Refs. p6|, [28|1). 



In the present case, certain similarities with a rotating black hole are 
found. The total flow velocity becomes equal to the velocity of light relative 
to the medium when 

When this equation is satisfied, 

^00 = 0. (52) 

For smaller values of r, ^qo becomes negative. The flow thus possesses an 
effective ergosphere lying inside the position defined by Eq. (51). 

The ergosphere will always exist for spiral flow. It extends outside the 
event horizon (if there is one), and also outside the apparent horizon. To 
verify the last statement, we calculate the total flow velocity: 

7V2 l + rV(fi2 + r2)" ^ ^ 

At the apparent horizon, we find from this equation that = l/n^ if 
i7 = and > l/n^ if f2 > 0. Thus, at the apparent horizon the total flow 
velocity is even for a very small vortex component larger than the velocity of 
light in the medium. This horizon is therefore located inside the ergosphere. 

We will call the outer limit of the ergosphere the static limit, in analogy 
to what is done in general relativity |^ . It will not be further discussed here 



(although we will briefly return to it in the final section). It is determined 
by the condition that the absolute flow velocity is equal to the velocity of 
light relative to the medium, i.e., Eq. (51). The static limit is marked by a 
vertical line in Fig. 4. 

5.3 Qualitative features of the light trajectories 

The horizontal lines marked by R2 — RA in Fig. 4 correspond to the horizontal 
lines i?2 — -R4 in Fig. 2. 
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There is no Region 1 in the "effective potential", meaning that there are 
no stable orbits for light. This is just analogous to the conditions outside a 
perfect sink, as discussed above. 

Region 2 is divided into Regions 2a and 2b in the figure. There is no 
qualitative difference between the two. They correspond to light entering 
the flow with an inverse impact parameter that at some point equals the 
"effective potential". Rays in Region 2 escape to infinity. 

Below the potential's bottom point, and above its top point, lie Regions 
3a and 3b. The former represents a photon with negative angular momentum 
captured by the flow singularity. The latter describes the same situation for 
a photon with positive angular momentum. It should be noted that Region 
3 may also describe a photon that is emitted somewhere outside the event 
horizon and escapes to infinity. 

Region 4 refers to photons emitted somewhere on the inside of the top- 
or bottom points of the " effective potential" . Note that the R4- lines cross 
the apparent horizon; there exists no condition that forbid photons to cross 
the apparent horizon from the inside. The event horizon, however, forbids 
any crossing so that the horizontal lines end there. 

Region 4 is divided into three subregions. Region R^a lies on the negative 
half of the the axis, whereas i?4b and i?4c he on the positive half. Two 
remarks should here be made: 

First, Region 4 contains only photons that are emitted in the inner re- 
gion of the flow and eventually become captured inside the event horizon. 
Although it remains adequate to keep the quantity 6 as a formal input pa- 
rameter, the interpretation of it as an 'impact parameter at inflnity' becomes 
lost. No photons in Region 4 will ever come from, or go to, inflnity. 

Second, the distinction between the lines R4i, and i?4c is not physically 
important. The first line corresponds to a photon with an inner turning 
point at b^^, while the second line has an inner turning point at bZ^. The 
two curves meet at the point where they are tangential to the apparent - 
horizon vertical fine. Together they enclose the inaccessible region of space, 
and there is no physical distinction between the two. 

Quite generally, the paths of light will be more complicated here than 
in Schwarzschild geometry. This is so because of the rotation. What can 
be read off from our plot of the "effective potential", is which trajectories 
are captured, and which can escape. In order to calculate the real shape of 
trajectories, the azimuthal equation of motion has to be solved. This task 
will not be undertaken in the present paper, but we will discuss the azimuthal 
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equation briefly in the next section. 



6 Summary, and final remarks 

The usefulness of picturing a non-uniformly moving medium as an effective 
gravitational field, endowed with a metric g'^", as demonstrated earlier for a 
pure vertex flow in Ref. turns out to hold true also for the present case 
with a spiral flow. The effective metric is defined in Eq. (6); its contravariant 
and covariant components are given in Eqs. (20) and (21). 

The equation for a photon in Fourier space, Eq. (4), is in terms of the 
effective metric picture expressible in the form of the Hamilton- Jacobi equa- 
tion, Eq. (5). At the turning point of a photon trajectory, dr/dX = 0. The 
turning point solutions of the Hamilton- Jacobi equation are given as solu- 
tions of Eq. (27), in which means the turning- point value of the inverse 
impact parameter b~^, in general defined as = E/L. There are in general 
two solutions of the equation, called and bZ^- The solutions are given by 
Eq. (29), and are identified as the "effective potential(s)" . They are shown 
in Fig. 2 for pure sink flow, in Fig. 3 for pure vortex flow, and in Fig. 4 for 
the general spiral case. "Effective potential" considerations turn out to be 
useful for the qualitative analysis of the various kinds of photon trajectories 
that can occur, in particular, for the general case of spiral flow. 

It is rather remarkable that the "effective potential" possesses properties 
that are qualitatively similar to those encountered for rotating black holes 
in general relativity. Thus, we can identify an event horizon, as defined by 
Eq. (45), for the flow if its sink component is strong enough in comparison to 
the vortex component to satisfy the condition kT"^ > Q^. A photon emitted 
on the inside of the event horizon can never escape to infinity. There exists 
also an apparent horizon, Eq. (49), lying outside the event horizon (if the 
latter exists). At the apparent horizon, the "effective potential" is cut off. 
Photons emitted beyond the apparent horizon, but not beyond the event 
horizon, can escape to infinity. Finally, there is an ergosphere, whose outer 
boundary is determined by the condition that the flow velocity is equal to 
the velocity of light relative to the medium, Eq. (51). In the ergosphere the 
time Killing vector C,{t) is spacelike (^oo < 0), so that no state of rest relative 
to infinity can be defined. 

We mention in passing that one interesting possible application of the 
pure vortex optical black hole is its relationship to the optical Aharonov- 
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Bohm (AB) effect This effect is in principle able to explore the long- 



ranging topological nature of a quantum vortex. The optical AB effect was 
analysed by Leonhardt and Piwnicki ||^, assuming a very high value of 
the refractive index. Because of the AB phase shift around a vortex one may 
say, strictly speaking, that we are considering angular momenta rather than 
impact parameters in the present paper. 

We ought finally to make clear what are the limitations of the present 
kind of theory. First, only nondispersive media are considered. Dispersion 
- sensitive effects such as the EIT - generated 'slow light' are thus not 
encompassed by the theory. 

Second, although the "effective potential" turns out to be useful for de- 
termining the turning points for photons, this concept does not allow us to 
determine the exact form of the spiral trajectories. If quantitative results 
are required, then the azimuthal equation of motion has to be solved. This 
equation can be constructed by considering the ^-component of Eq. (18): 



and solve the equation 



^ = -g'^'K = -g'^X + g^^'K + g^%- (55) 

If this is done, the actual trajectory can be plotted. A 0-diagram would show 
special features at the point where the azimuthal flow component becomes of 
the same magnitude as c/n. We have made some calculations in this direc- 
tion, but the formalism becomes somewhat unwieldy and is not reproduced 
here. 



6.1 On the physical meaning of the horizon 

It is of interest to carry out the discussion of the physical meaning of the 
horizon a bit further. In general relativity, a horizon has interesting quantum 
properties that have their roots in the behaviour of classical waves. This 
is related to the phase of a wave developing a logarithmic singularity at 
the horizon. The well known Hawking radiation stems from a Bogoliubov 



transformation relating the "in" vacuum to the "out" vacuum One may 
wonder: does our physical system allow Hawking radiation to occur? 
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To investigate this point further, let us first recall some facts about the 
behaviour of a scalar field $ near the horizon r = 2M of a Schwarzschild 
black hole (for simplicity we now put G = c = 1). As shown by Jensen 
and Candelas ||2^, near the horizon the solution of the scalar field equation 



□$ = 0, regular at infinity, varies with r as 

/ r \ —2iMu) 

$ oc ( 1 ) . (56) 

V2M / ^ ' 

The wavelengths become shortened near the horizon; the waves pile up. The 
effect is demonstrated graphically, for instance in Ref. (dealing with the 
somewhat more general case of a Schwarzschild - de Sitter scalar field). It is 
worth noticing here that the effect is quite analogous to that taking place near 
the horizon in a Rindler frame. In that case, the field mode corresponding 
to Eq. (56) is found to vary as the modified Bessel function Ki^ of imaginary 
order [||] 

^(xKUkx). (57) 

Here k = |k^| is the magnitude of the transverse wave vector k^, and x is 
the longitudinal distance from the horizon. Since Ki^{kx) = ^T{iu){^kx)~''''^ 
for small arguments, it follows that $ oc x"*"^, thus in essential agreement 
with Eq. (56). 

After these preliminaries we return to our fiow field. To avoid complicat- 
ing the formalism too much, we imagine a scalar field $ that is azimuthally 
symmetric. In the effective metric picture the field equation is 



d, = 0, (58) 

where the indices /i and u run over and r. From Eqs. (21) we calculate 
the metric 
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1--T 1 + 



^2 V J.2 



,2 + + 



We will investigate whether there are solutions of Eq. (59) that are singular 
at the horizon. As "horizon" it is natural here first to consider the static 
limit, r = r stat.limit, since ^oo = there. According to Eqs. (51) or (21), 



^stat.limit 



y/K{n^ + T^). (60) 
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Consider now positions r just outside the static limit. As the metric ap- 
proaches a constant value in this limit, g — > — + F^), we can omit 
the influence from this factor altogether when looking for singular solutions, 
and consider instead the reduced equation 

Cg'^-'d,^) = 0. (61) 

Assuming a time factor exp{iujt) we can write this equation approximatively 

as 

+ ^ ^ tu^$ \ 1 H — dr^ 



/ 7^ V 



0. (62) 



This equation actually has no solution singular at the static limit. Such a 
solution occurs, however, at the apparent horizon, r = Vapp = \fK,\T\. The 
first term in Eq. (62) can be neglected to leading order, and the scalar mode 
vanishing at r ^ cxd becomes, close to the apparent horizon. 



<|.cx(l-^ . (63) 



This field, in fact, has the same qualitative behaviour as that outside a 
Schwarzschild black hole, Eq. (56). Physically, this singular behaviour can 
be understood from the fact that the apparent horizon is the point of no 
return for light coming in from infinity. 

Does this imply that there really occurs Hawking radiation from the op- 
tical black hole? The answer is in our opinion no, the reason being that 
the metric cjf^i, is an artificial metric only. Our fluid system is physically a 
special-relativistic system. There exists one symmetry group in the present 
case, namely the Poincare group, according to which the vector d/dt is a 
Killing vector of spacetime (cf. page 45 in |21|]). The condition for Hawking 
radiation to occur, is that there are physically inequivalent "in" and "out" 
vacua. That is the case in curved spacetime, where the Poincare group is no 
longer a symmetry group. (This point is however not quite trivial, and it is 
only fair to mention that Hawking radiation has been claimed to occur in 
recent discussions on acoustic black holes.) 
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